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The  benefits  of  feedback  achievable  in  nonminimum-phase  (nmp) 
systems  are  severaly  limited.  However,  if  there  are  two  parallel  nmp 
branches  whose  outputs  can  be  measured,  then  it  may  be  possible  to 
eliminate  or  reduce  the  nmp  property  with  respect  to  some  parts  of  the 
system,  even  if  each  of  the  two  branches  has  uncertain  parameters.  It 
is  always  possible  to  do  so  for  a finite,  nonzero  range  of  uncertainties 


) 


for  any  complexity  of  pole-zero  patterns  of  the  branches.  A design 
philosophy  and  methodology  is  presented  for  the  general  problem,  permitt- 


ing one  to  approach  an  optimum  solution,  under  various  constraints. 
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SYMBOLS  AND  ABBREVIATIONS 


B(ui) 


mp 

mps 


NL 

nmp 

rhp 


boundary  of  T (u>)  (2.1) 

number  of  encirclements  of  origin  by  tfr  (2.) 
minimum-phase  ( 1 . ) 
mp  and  stable  (2.1) 

boundary  of  quadrant  1 (2.)  , map  of  F on  . 

Nyquist  locus  — map  of  a function  on  (2.) 

nonminimum-phase  ( 1 . ) 
right  half  plane  (1.) 


T (w)  {P(jw)}  (2.1) 

0p  Arg  P(j«a) 
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bunding  of  uncertain  nonminimum-phase 

PLANTS  FOR  ELIMINATION  OR  REDUCTION  OF 
NONMINIMUM-PHASE  PROPERTY 


1.  INTRODUCTION 

Consider  a feedback  loop  containing  a constrained  part  (denoted  as  the 

Plant)  with  transfer  function  P(s)  , and  compensation  G(s)  to  be  chosen 

T a 3t /t 

by  the  designer.  Let  Sp(s)  £ — ^ be  the  sensitivity  of  any  system  input- 
output  function  T(s)  to  the  plant  P(s)  . It  is  known  (Horowitz  1963)  that 


sp  - Tfir  (1-t}  ' where  T0 


(Tjp^Q  is  the  'leakage  transmission*.  Hence 


Lp (s)  which  is  the  loop  transmission  for  reference  P " (Horowitz  1963)  is 

of  great  importance  in  feedback  design,  especially  in  the  large  class  where 

T 1 

the  plant  output  is  the  system  output,  for  which  Tq  ■ 0 , giving  Sp  * 1+L^  * 

If  P has  no  zeros  in  the  right  half-plane  (rhp) , i.e.,  is  minimum- 

phase  (mp) , then  theoretically,  at  least,  | Lp( jw) | can  be  made  as  large  as 

desired  over  any  finite  bandwidth  (Horowitz  1975) . However,  if  P is  nonminimum- 

phase  (nmp)  then  Lp  is  severely  constrained  (Horowitz  and  Sidi  1978) . For 

example,  suppose  P - Pm(s-a)  , a > 0 , Pm  mp  and  Lp(j<a>)  is  designed  in 

the  usual  manner  with  | Lp  ( ju>)  | > 1 for  u < o>c  , <1  and  monotonically 

decreasing  for  u>  > u . Then  (u  ) *•  0.5a  , i.e.,  the  available  feedback 

c c max 

loop  bandwidth  and  with  it  the  benefits  of  feedback  are  severely  curtailed. 

It  is  therefore  very  important  to  eliminate,  if  possible,  or  at  least  alleviate  the 
nmp  property  in  a plant.  It  is  shorn  in  this  paper  that  to  some  extent  and 
for  a certain  class,  this  can  be  done.  Note  that  rhp  poles  of  a plant  do  not, 
by  themselves,  limit  the  feedback  bandwidth  achievable  (Horowitz  1963). 


2- 


' 

I 


I 
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1.1  A Class  of  Plant*  for  which  nmp  Alleviation  is  possible. 

In  Figure  1,  the  constrained  part  (darker  lines)  has  a single  input  x 
and  two  outputs  , y2  independently  neasurable  and  compensa table  by 
Ga  , . The  transfer  functions  Yj/X  , y2/X  have  Pj  in  comnon  (usually 

the  poles) . The  transfer  function 


Z(s) 

X(s) 


P.  (P  G 

l a a 


*W 


(1) 


Suppose  P^  , Pfa  are  nap  and  with  uncertain  parameters.  Is  it  possible  to 

find  Ga  , G^  such  that  (paGa  ♦ Bp  over  tlie  ®n^ire  range  of  these 

uncertain  parameters?  If  so,  then  the  loop  transmission  for  reference  P^  has 

no  rhp  zeros  due  to  P#  , . If  is  also  mp,  then  the  system  mav  be 

designed  to  have  arbitrarily  small  sensitivity  to  PJ  . In  any  case,  the 
T 

limitations  on  S due  to  the  nmp  P_  , p.  are  at  least  eliminated.  However, 

a d 

it  is  important  to  note  that  L#  , the  'loop  transmission  for  reference  P^  ' , 
is  obtained  by  cutting  the  loop  at  aa'  and  so  does  have  the  rhp  zeros  of 
Pa  , and  similarly  1^  must  have  those  of  Pb  . The  incremental  linear 
model  of  the  longitudinal  axis  of  an  aircraft  is  an  example  of  Figure  1,  with 
x the  elevator  control  surface,  yx  pitch  rate  and  y2  normal  acceleration 
at  pilot  station.  Under  certain  conditions,  P and  P^  are  both  nmp  and 
with  highly  uncertain  parameters  (Edwards,  Rediess  and  Taylor  1970).  Other 
examples  are  common  in  process  control,  with  time  delays  often  constituting 
the  nmp  paths.  These  may  be  well  modelled  by  all-pass  networks  over  any 
desired  frequency  range. 

There  is  no  single  or  set  of  formulae  available  for  this  problem. 

Rather,  this  paper  presents  a design  philosophy  and  methodology  usable  for  any 
specific  case.  It  is  developed  here  by  proceeding  step  by  step  from  simple 


to  more  complex  problem  classes 
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2.  THE  SINGLE  POLE-ZERO  PAIR 

In  Figure  1,  the  zeros  of  concern  of  (PaGa  + Pb<l  ) are  those  of 

f G Pv.  G Pk1 

-+-b+p5*(  F » , P « rM  (2) 

<■  °b  Pa  °b  a J 

F is  to  be  chosen  so  that  ♦ has  no  rhp  zeros  over  the  range  of  uncertainty 
of  P . Considerable  insight  is  obtained  by  studying  the  simple  class 
K k(s-z) 

P * -7 — r- r - £ K.JcP.  (s)  , K,k,z  > 0 , X > 1 , with  k € (l,k_],  k,  > 1 

(S-AZ)  111  i 

denoting  the  uncertainty  only  in  k in  the  meantime.  Hereafter  ■ 1 is 

used,  so  that  in  practice  the  F that  emerges  must  be  divided  by  the  actual 

Kx  . Note  from  (2)  that  if  F is  a solution  for  an  uncertain  set  {p}  , then 

1/F  is  a solution  for  the  set  tl/P}  • Also,  F need  not  have  an  excess  of 

poles  over  zeros.  Let  be  the  boundary  of  half  the  forbidden  region  — the 

upper  rhp  (first  quadrant).  By  the  Nyquist  locus  (NL)  of  H(s)  denoted  by 

H^  , is  meant  the  mapping  of  H(s)  as  s moves  clockwise  on  from  0 

to  jR  and  then  s * Re^0  , 0 from  s/2  to  zero  , R -*•  » . P^  is 

sketched  in  Figure  2a  at  k ■ 1 and  k ■ k2  < X . It  is  easy  to  see  that  if 

ko 

-F  is  a constant  € (— » 1)  , then  V k 6 (l,k.)  , of  (2)  executes 

x 2 % 

half  a negative  encirclement  of  the  origin  (denoted  by  ■ -0.5)  . Invoking 

the  Nyquist  criteria:  N^  - No.  of  encirclements  of  origin  - 2£^  - -Np 

where  N - no.  of  rhp  zeros  of  ♦ , N_  * no.  of  rhp  poles  of  ♦ - 1 here, 

giving  N ■ 0 . 
z 

However,  a dynamic  F(s)  is  needed  if  k2  > X (Figure  2b).  The 
restriction  In  F regular  in  rhp  is  made  in  the  meantime,  for  which  the  Bode 
integrals  (Bode  1945,  Horowitz  1963)  are  very  helpful  in  deriving  sup 
optimum  F , denoted  by  Fgup  ~ since  a strictly  optimum  F does  not  exist 
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but  Fgup  can  be  approached  as  closely  as  desired.  It  is  next  shown  that 

2 

Fsup  “ (s+Az)/(s+r)  , giving  sup  k2  * A » with  the  root  loci  of 
1 + JcP * P ■ JcP^  , shown  in  Figure  4 . 

2.1  Derivation  of  F 

....  .<■  — sup 

Let  T(u>)  (template  of  P(ju))  ■ {P(Jw)}  over  the  range  of  uncertainty 
of  P , and  B(u>)  ■ Boundary  T(w)  . If  the  uncertainty  is  only  in  k then 
T(u>)  = B(w)  , sene  of  which  are  shown  in  Figures  2b,  3a,  b the  numbers  and 
letters  in  order  of  increasing  u values.  The  NL  of  a satisfactory  -F  is 
shown  in  Figure  2b  and  in  Figure  3b  on  the  Nichols  chart.  (Henceforth  the 
subscript  % 1*  omitted  in  the  figures.)  F(s)+P (s)  is  the  complex  number 
(or  vector),  originating  at  -F(s)  and  terminating  at  P(s)  . By  considering 
the  motion  of  P relative  to  F , i.e.,  taking  -F(s)  as  the  origin  of 

(F+P)  at  each  s € % , it  is  seen  that  P encircles  -F  one  half  times 
counterclockwise,  ■ - .5  , V k € [l,)^]  • Mote  that  Arg(-F)  > Arg  P 

in  Figure  2b,  except  at  w • 0 , • and  that  a -F  similar  to  above,  but 
with  Arg(-F)  < Arg  P would  give  £^m  *5  * which  is  not  acceptable.  It 
is  seen  that  the  following  are  necessary  and  sufficient  conditions  for  F+P 
to  be  mp  V k € U, k2)  . 

A.  (a)  -F(0)  > k2P(0)  . (b)  -F(-)  < P(-) 

B.  -F(s)  * T(s)  V s € (3) 

C.  ^ ■ - .5  , with  the  following  geometric  interpretation: 

In  Figure  3b  the  B(w)  are  vertical  lines  of  equal  length.  Each  point  on 
B(u)  corresponds  to  a specific  k value  , Bk(w)  . At  any  fixed  k , 

let  n be  the  timber  of  u values  at  which  |F(ju)|  - | Bk(w)|  , 
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I 


Arg  -F(jw)  < Arg  B^tw)  . Then  n must  be  even  in  order  that  * " *5  , 

as  in  Figures  3a,  b.  In  other  words,  -F(jw)  in  moving  from  above  B(0) 
to  below  B (oo)  , as  per  (3Aa,b),  may  temporarily  be  on  the  left  of  B(w) 

(Arg  -F  < Arg  P)  in  Figures  (3a,b) , but  in  the  net,  must  execute  the  passage 
on  the  right  of  B(o>)  where  Arg  -F  > Arg  P . 


From  (3Aa,b),  it  is  necessary  that 


the  maximum  tolerable  k 
maximize  |f(0)/F(«)|  , 
(denoted  by  mps)  F(s) 


in 


F(0) 

F(») 


(recall  k^  - 
subject  to  (3B,C) 
(Horowitz  1963) 

=°  6 _(u)dd) 

— , 0 . 

u -F 


F(0)  P(0)  _ ^2 

F(-)  * 2 P(»)  X ' 


with 


1) . Hence  it  is  desirable  to 
. However,  for  mp  and  stable 


(u)  Arg  [-F(jw)] 


k 


2 


(4) 


so  it  is  desirable  to  minimize  0_p(u>)  at  each  m . 

If  the  above  defined  n « 0 for  all  s € ft*,  then  from  (3B) 

0_F(«o)  < Arg  B(u>)  - Arg  P(ju>)  & 0p(o>)  , V u € 10,-)  , giving 

inf  0 _(w)  - 0_(<i»)  . Since  F is  mps,  it  is  determined  (up  to  an  arbitrary 
•F  P 

constant  factor)  by  0p  , giving  the  supremum  F(s)  , denoted  by 


_ , . -(s+Xz) 

Fsup  " “1/P!(*)  " K(s+z) 


2 2 2 2 2 

K is  obtained  from:  Numerator  (F  +P)  - (s  -X  z ) -kX(s  -z  ) - 

sup 

- s2(l-kK)+z2(kK-X2)  . The  coefficients  must  have  the  same  sign,  giving 
kjK  — K — 1 , k2K  - X2  - k2  . The  root  loci  of  1+kPj/F^  ar®  sketched  in 
Figure  4.  *gUp  is  obviously  unique. 

2.1.1  Case  n f 0 . Obviously  from  (4)  it  appears  desirable  to  let  0_p  < 0p 
for  some  intervals,  e.g.,  in  (0,o»x)  in  Figures  3a, b with  0_p  “ 0p  fgr 


4 
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u € (.  ,-)  . The  value  of  (4)  would  be  presumably  thereby  Increased.  The 
trouble  is  that  (3C)  dictates  |Kj<*x)|  > Jt^ I ^i^ J**x>  I with  |f(»)|  < |p(*)|  * 1 

I 1 P(j“x) 

as  before,  giving  k2  < j - ^ f . But  increases  with 

io  , as  does  | F («) /F  ( ju  ) | . as  shown  below.  Thus  (Horowitz  1963) 
x x min 


i p(>*>  ( r 2 

*n  I “FI=0“|  " / w(w)  0-F(“)d“  ' *•  - 7 -f- 


.2  2 ' 
(u>  - u ) 

X 


with  ji(w)  >0  for  w < ux  , <0  for  w > o»x  . Hence,  to  maximize  (5)  , 

maximize  6 in  u < w , minimize  it  in  u > u . However,  the  set 

“F  X X 

{F(ju>)}  is  constrained  to  have  6 < 0_  for  u < u>  , 0 > 0 in 

*F  r X *F  r 

(ii  > «i»x  , so  let  sup  0_p,  “ 0p  V W . Thus  (5)  decreases  vs  o>x  , and 
(“x)opt  " 0 ' 9ivin9  the  previous  Fsup  . 

with  General  Uncertain  Poles  and  Zeros. 

In  P - k(s-z)/(s-p)  , k € tl,k2J  , suppose  z and  p have  independent 

uncertainty.  Then  a typical  T(w)  is  shorn  in  Figure  3c  for  z < p . 

max  mm 

B (u)^  the  boundary  of  T(«) ^ has  vertical  lines  due  to  the  independence  of 

the  uncertainties.  It  should  be  clear  from  2.1,  2.2  that  the  worst  case 

determines  max  k2  to  be  (- 1 . This  is  also  seen  by  considering  only 

the  P . , z case,  for  which  k_  has  the  above  value.  The  resulting 

min  max  2 max 

design  also  obviously  satisfies  all  other  p , z smaller  k2  combinations. 

If  P . < z and  P > z . , then  nap  elimination  is  impossible  for  any 

min  max  max  min 


finite  k uncertainty,  no  matter  how  small. 
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Consider  P as  in  above,  but  uncertain  z,  p,  k related  to  some 
extent-,  giving  more  general  T(u)  . The  previous  arguments  leading  to 
-F  lying  on  B(w)  V u>  € still  apply.  The  proof  of  existence 

Sup 


and  uniqueness  of  such  an  FSUp  is  basically  the 


as  in  (HOrowitz  and 


Sidi  1978).  However,  it  is  more  difficult  now  to  find  F , because 

sup 

unlike  the  above  Arg  Fgup(w)  is  not  apriori  known.  It  is  only  known  that 
there  exists  a unique  mps  analytic  function  which  lies  on  B(u>)  at  each  o>  . 

( S“Z  ) 

Finally,  consider  P * k y P2(s)  , P2  having  only  thp  poles  and 

(s+a^) 

zeros  with  uncertainty,  e.g.,  P2(s)  ■ Ts+a  ) (s+a  ) ' ai  € ^ail'ai2*  * In 
order  that  the  new  T(w)  and  B(u)  be  qualitatively  similar  to  the  previous, 

Fl<s+ao> 

we  could  let  F ■ ; — d ■/  ■ — r-  , a fixed  in  [a. ,a_]  etc.  and  then  the 

o o 

problem  is  one  of  zeros  of 

<_  ....  (s+b)  (s+c„) 

F + v SlZlL  .1*+*).  2-  2_  A f +p’ 

1 (s-p)  <s+aQ)  (s+b)  (s+c)  * *1+F  * 

Again,  the  reasoning  in  2.1,  2.2  applies  leading  to  the  conclusion  that  Fsup 
lies  on  B (ui)  at  each  u>  € [0,»)  , where  B (w)  is  that  of  P'  over  its 

range  of  uncertain  parameters.  In  the  balance  of  the  paper,  only  uncertainty 
in  k is  considered,  with  generalizations  applicable  in  the  same  manner  as 
the  above. 

2.4  More  Conservative  Boundaries. 


It  is  obvious  that  the  above  design  philosophy  applies  to  other  kinds  of 
bounds  on  the  zeros  of  F+P  ■ i>  . A reasonable  alternative  is  that  the  zeros 
of  <l>  must  have  a real  part  t - o , a > 0 . The  line  Re  s • - o replaces 
the  imaginary  axis,  so  let  s - w-o  , P(s)  - k giving 

* and  resulting  new  k ■ X2 < X2  and  new 

new  o+z  ’ sup  new 
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w+(o+Xz)  s+2p+Xz 


sup  w+(o+z) 


s+2o+z 


(6) 


Another  possibility  is  that  "damping  factor  of  the  zeros  of  >,  > 0 . 

The  boundary  is  then  a straight  line  in  the  second  quadrant  through  the 

origin,  at  sin  * £ with  the  vertical.  This  is  the  new  % but  the 
o 

procedure  and  results  are  as  before,  i.e.,  one  finds  B(s)  , s , etc., 

and  develops  the  integrals  analogous  to  (4,5)  in  order  to  derive  F 

sup. 

2.5  Reduction  of  Nonminimum-Phase  Property. 

2 

In  2.1,  2.2  if  kj  > ) , then  for  some  k , $ of  (2)  must  be  nmp. 

However,  one  may  design  so  that  the  'worst*  nmp  < that  in  an  uncompensated 
design.  The  concept  of  less  or  more  of  the  nmp  property  is  evident  by  noting 
that  rhp  zeros  contribute  excess  phase  lag  (Horowitz  1963) . The  further  the 
rhp  zero  from  the  origin,  the  less  this  excess  lag  vs  u , so  if  the  nmp 
property  cannot  be  eliminated,  it  is  desirable  to  at  least  guarantee  that  the 
rhp  zeros  lie  as  far  away  as  possible  from  the  origin. 

In  this  new  problem,  \ - WjW2...W5  , semicircle  of  radius  R in 
Figure  5,  because  the  forbidden  region  for  the  zeros  of  i>  is  its  interior. 

But  the  reasoning  and  technique  is  precisely  the  same  as  before.  Again,  the 
conclusion  is  that  FQpfc  (not  sup  because  the  boundary  of  fy  is  assumed 
permitted)  lies  on  B(s)  for  each  s € % , and  this  corresponds  to  each 
s €^on  a root  locus  of  FQpt  + for  k € (l,k2)  . 

The  derivation  of  Fopt  is  not  difficult  for  the  case  of  k uncertainty 
only.  Arg  Pj/P  t must  have  a constant  value  on  (zero  or  180",  the 

former  here) . The  inverse  image  with  respect  to  v»2  W3  , of  the  zero-pole 
pair  at  z , Xz  is  a pole-zero  pair  at  R /Xz  , r/z  , i.e.,  this  combination 


T 
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! 

t 

I 


t 


of  2 pole-zero  pairs  has  a constant;  angle  on  W2  W3  • The  inverse  image 
of  these  two  pairs  with  respect  to  is  their  ihp  image,  shown  in 

Figure  5,  with  the  resulting  root  loci.  The  new  )l2max  is  now 


2 2 2 2 

\'(K-S)  > x2 

,2  .2  2.2 

(R  -X  z ) 


e.g.,  if  z - 2,  X - 4 , R - 10  , then 


(7) 


ne”  k2ux  64 

jjTj— --T.  "ere  Fopt  h»s 


rhp  pole  and  zero,  apparently  violating  the  constraint  that  F is  mp  and 
stable.  This  is  because  the  actual  constraint  used  is  that  InF  be  regular 
in  the  forbidden  region,  which  is  now  the  inferior  of  W^. . .w5  . 

Of  course,  one  may  define  different  shapes  for  ?£/  , e.g.,  in  a 
specific  problem,  a highly  underdamped  rhp  zero  pair  may  be  more  tolerable 
closer  to  the  origin  than  a more  damped  pair  further  from  the  origin 
(Horowitz  1963) . An  ellipse  with  horizontal  major  axis  might  then  be  more 
suitable  for  . The  conclusion  that  *opt  lies  on  B( 


for 

s € 

Fopt  ' 

or 

good  approximation  to  it. 


- 
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3.  MORE  COMPLEX  P(s) 

While  the  results  in  Section  2 are  of  interest  in  themselves,  they  are 

more  important  in  suggesting  an  approach  usable  for  more  complex  P . An 

important  conclusion  was  that  -F  (u>)  € B(w)  , V u € [0,®)  which  coincided 

sup 

with  0_p  * 0p  . The  more  complex  problems  next  treated  show  that  in  most 
cases  0 _ ■ 0_  doe3  not  coincide  with  -F  € B(u>)  V u . Each  case 

-F  P SUp 

considered  introduces  a new  property. 

3.1  Case  1 - High  Frequency  Modification. 

Consider  P - A k P (s)  , k € [l,k_]  with  p,E,H  > 0 , and 

s -Es+H 

such  that  0p  €(-.5ir,0)  in  (0,u*)  , 6(0,. 5x)  in  ((d#,«)  and 

|P(jw) | > P(0)  , as  in  Figures  6a,  b.  If  k2  < 0Q/0A  in  Figure  6a,  then 
a constant  value  for  -F  with  k^tO)  < -F  < 0Q  » P^jw*)  is  satisfactory. 
But  for  larger  k2  a dynamic  F(s)  is  needed,  such  that  -F  encircles 
{B(u) } once  negatively.  Given  -F(0)  > K2P;l (0)  , this  requires  the 
existence  of  an  u>x  at  which  |F(jwx)|  < IP^jw^l  , e_F(“x)  ■ eP<“>x)  * 

«md  -F(»)  > k2P1  (•)  ■ 0 . (It  is  easily  seen  that  -F(0)  < Pj^  (0)  is 
inconsistent  with  mps  F .)  Also  in  Figure  6b,  -F  must  in  the  net  execute 
the  passage  from  above  B(0)  to  below  B(ux)  on  the  right  of  B(«o) 

(0_p  > 6p)  , and  the  passage  from  below  B(«x)  to  above  B(«»)  in  the  net 
on  the  left  of  B(o»)  with  0_p  < 0p  (cf  2.1  after  (3)  and  Figures  3a,  b) . 

An  acceptable  -F  is  shown  in  Figure  6b  with  OKI  , OK2  two  among  many 
possibilities  for  w > u>x  . 
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The  condition  on  F at  “ is  easy  to  achieve , so  k2  maximization 


involves  maximization  of 

2w2 


|F(0)/F(jo>x)  | - — - / w(o))e_pda)  ,*y(fa))- ^ 


(8) 


fa)  (fa)  -fa)^) 


A 

o 

for  o>  < u>x  , >0  for 

fa)  > fa>x  . Hence, 

it  is  desirable 

to  minimize 

CD 

1 

in  (0  < wx  and  maximize 

it  in  fa)  > fa>x  . 

However,  it  was 

noted  above 

that 

in  the  net  0 _ > 0_  in 
-F  P 

(0,0)^  . 

One  might  suggest  use  of  0_p  < 0p  for  some  intervals,  i.e.  n f 0 
as  in  2.1.1,  but  a technique  similar  to  that  of  Equation  (5)  can  be  used  to 
prove  it  worse  here  also.  Thus  it  would  require  existence  of  an  < u>x  at 
which  0_F  ■ 0p  and  |f|  > k2)PjJ  . Maximization  of 


2 2 

_ « fa>0  (w  -w  )dfa> 

tn  |r(j.0)/F(J.  )|  --/  - » ° 

0 (fa)  -0)  ) (fa)  -fa)  ) 


derivable  from  (8)  by  considering  tn 


F(  jw  ) 
x 

F(0) 


(9) 


and  In 


F(jwo) 

F(0) 


leads. 


as  in  2.1.1,  to  u - 0 . 

o 

The  above  suggests  trying  9_p  - 0p  , i.e.,  -F  - -F^  - — , 

qualitatively  shown  in  Figure  6a.  It  is  interesting  to  relate  this  to  the 
root  loci  of  1+P/F^  shown  qualitatively  in  Figure  6c.  The  numbers  E - 5, 


F ■ 6 , p * 1 are  taken  for  illustration.  The  double  root  at 

- j3.9  corresponds  to  M (on  P . ) - note  that  -F. 

A min  1 

lies  on  B(u)  on  both  sides  of  M in  Figure  6a  and  to  the  lowest  k for 


which  the  root  loci  are  all  on  the  ju  axis  in  Figure  6c.  The  double  root 


-12- 


at  the  origin  corresponds  to  V 
-F  when  the  locus  for  w < 0 


(at  k2)  corresponds  to  point 
B(w)  only  for  u < u . 


in  Figure  6a  — with  a similar  property  of 
is  included.  The  single  root  at  j4.8  =* 

Y in  Figure  6a  — at  which  -F  lies  on 


Since  -F^  lies  on  B(u>)  for  u < u>x  , no  improvement  (decrease  of  9_p) 

is  possible  there.  But  it  is  possible  to  improve  matters  by  increasing  e F 

for  u > m because  -F,  does  not  lie  on  the  B(u)  for  u > u . Any  such 
y l y 

increase  will  however  decrease  |F/F(0)|  V 0 < u < for  some  A > 1 

(Horowitz  1963) , creating  new  id^  , and  larger  . A simple  way  to 

(s2+Es+H)  (s2-hd2) 

increase  0_p  for  <»  > > <oy  is  to  let  -F2  - i(s+'p) ~ * which 

leaves  0_p  unaltered  in  (O.u^)  . It  is  relatively  simple  to  find 

<i>z  ( - 10.9)  from  Num.  (F2+P)  by  requiring  that  at  kj  « 1 , it  has  a factor 
2 2 2 

(s  +01^)  i Ul  - 2.1  replaces  in  Figures  6b, c.  At  k2  it  should  have 

2 2 2 

a factor  s (s  +u>2)  giving  u>2  ■ 7.24  replacing  . (Compare  Figures  6a, c 
at  w ■ 3.9,  7.24,  9.8,  10.9).  ~F2  9V  *8  s*K>wn  *-n  Figure  6a  and  permits  a 

slightly  larger  k2  - 1.797  in  place  of  the  1.731  for  Fx  . Note  that  -F2  lies 
on  B(u)  for  w 6 (0,9.8)  > (0,4.8)  for  which  -F^  lay  on  B(u)  . 

One  may  increase  k2  by  repeating  the  above  with  the  factors 


(s2+a 


m+1 


) 


Ai  < Bi  < *i+1  . Such  increases  are  always  possible  if  3 any  interval  in 
which  -F  does  not  lie  on  B(«)  . However,  the  increase  obtained  is  found 
to  be  exceedingly  small  and  does  not  justify  the  design  effort. 
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3.2  Case  Two  - Low  Frequency  Modification. 

2 2 

Let  P - k(s  -As+B)/(s  -Cs+D)  6 k Pj  , k € {l,k2>  with  parameters  such 

that  P is  as  shown  in  Figure  7a,  i.e.,  .5*  > 0p  > 0 in  (0,u^)  , 

- .5ir  < 6p  < 0 in  (o^,*)  ; | P ( jtu^)  | < |P(»)  | < P(0)  . Again,  a constant 

-F  can  handle  uncertainty  k2  < |P^  ( | 1 , if  k^P^jw^)!  < -F  < 1 . 

For  larger  k2  a dynamic  F(s)  is  needed  satisfying  the  following  conditions 
-F(0)  < P1(0)  ; F(«)  < Px (~) 

3#x  3 0p  * e_F  but  -F(jo)x)  > k2  P(ja>x)  | (1°) 

F(s)  t k P-^s)  , s € % and  £ * - 1 . 

From  (10)  it  is  desirable  to  maximize  |F(j<i>x)  |/|f(0)  | , |F(  jwx)/F(<») 

The  effect  of  e on  the  former  is  given  by  (8).  Its  effect  on  the  latter 
is  (Horowitz  1963) 


in 


F(w  ) 
F(«) 


■ —<i)  / 0 _(u)A(«)du  , A(w) 

* * 0 ~F 


u 


2 2 
u -u 

X 


(11) 


Thus,  from  the  signs  of  the  weighting  functions  p(u)  , A(u)  in  (8.11),  it 
is  desirable  to  minimize  0 € (0 ,w  ) and  maximise  it  in  («,«•)  subject 

“F  X X 

of  course  to  (10).  However,  from  previous  experience  one  might  guess,  and 

it  is  in  fact  so,  that  in  these  intervals  inf  and  sup  0_F  are  equal  to 

0p  , so  it  is  reasonable  to  try  (s  +Cs+D)/(s  +A8+B)  4 -Pj  — Figure  7a, 

OB 

giving  k2  * — . It  is  seen  that  u>x  of  (10)  • and  tnat  indeed  6p 
constrains  e p but  only  in  (u^,«»)  because  -F^  lies  on  B(u)  only  in 
(o^,-)  . Note  that  n«*0  is  assumed.  Proof  that  is  inferior,  is 


similar  to  that  previously  given  in  2.1.1,  3.1 
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Accordingly,  0_p  may  be  Increased  in  (0,<u£)  (because  the  change 

changes  u^)  , but  in  such  a manner  that  0_p  ■ 0p  in  . A simple 

s2 

way  to  do  this  is  to  let  F2  * F3  — ^ ~ • giving  0_p  *»  0p  for  ui  > u>2 

1 (s  +u>2) 

2 

but  0_F  ■ Op  for  u < u>2  • The  value  of  u>2  is  found  by  considering  the 
numerator  of 

. . P , Jc(s2-As+B)  (s2+As+B)  (s2+w2) 

1 + p * 1 ' 2 — 2 2 

s (s  -Cs+D) (s  +Cs+D) 

whose  (qualitative)  root  loci  are  shown  in  Figure  7b,  for  the  case 

P ■ I-8-  2>.  . One  chooses  a»_  so  that  at  k • 1 , the  numerator  has  two 

(S" |S“3/  2 

2 2 2 

zeros  at  «•  and  a factor  (s  +w  ) . Note  that  -F,  lies  on  B(u)  only  for 

o l 

u > Wjj  ■ .64  in  this  example,  while  -F2  lies  on  B(w)  for  u > u>a  ■ .J  . 

Note  the  relations  between  j>  , (i)  , (i)  and  the  k values  in  Figures  7a,  b 

o a c 

cf  similar  comparison  of  Figures  6a,  6c.  The  difference  in  k2  is  however, 
not  large  — 1.41  ■ OB'/OA  instead  of  1.39  - OB/OA  in  Figure  7a.  Increased 
iaproveamnt  is,  of  course,  possible  because  -F2  does  not  lie  on  B(u)  C (0,.3), 
but  it  is  already  seen  that  the  improvement  will  be  very  slight  — for  example 
use  of  (s2^)2/.4  gives  V ^ 1.412  with  -F3  on  B(u>)  € (0,.138)  , (.287,-)  > 

(.3,»)  for  -F2  . 

3.3  Case  Three  - High  Frequency  Modification. 

Zn  this  case  P has  a pair  each  of  rhp  poles  and  zeros  3 0 < 0p  < .5* 

in  (0,wo)  , - .St  < 0p  < 0 in  (u>0,«»)  # lp(3<*»0)|  * p(0)  * p(*>  * e.g., 

2 2 
p ■ k , k € tl,k9]  in  Figure  8a.  using  -F.  - — --y*3— 4--  gives 

(s  -3s+4)  A (sSs+2) 

K • 4 kj  » 2.88  and  -F  on  B(w)  for  w € (0,2.08)  . It  is  therefore 
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> 


! 

I 

1 


i 


possible  to  improve  0_F  in  the  high-frequency  range.  Here,  the  constraints 
on  a realistic  -F  are:  F(0)  < P(0)  and  at  some  , 0_p  ” 0p  with 

| F ( | > X2 1 P ( jtu^)  | , so  maximization  of  |F( jw^/FfO)  | is  desirable. 

From  (8)  this  is  achieved  by  maximizing  0_p  in  (0,u>^)  and  minimizing  it  in 
. But  such  a -F  will  be  on  B(w)  in  (0,w^)  , so  the  iitprovement 
can  be  made  only  in  (o^,-)  . This  is  conveniently  done  by  means  of  poles  at 
± jw,  * which  leaves  0 * 0 in  (0,«.)  . The  resulting  root  loci  of 

1 + P/F2  are  then  shown  qualitatively  in  Figure  8b,  for  the  above  numerical 
example.  To  optimize,  is  chosen  so  that  k^  ■ k^2  ■ k^  * 1 . The  result 

is  k2  = 3.23  > previous  2.88  with  ■ 3.987  , and  -F2  lies  on  B(w)  for 
u € (0,3.76)  , cf  previous  (0,2.08)  . In  this  case  use  of  two  coincident 
complex  pole  pairs  in  place  of  the  single  pair  leads  to  small  increase  of 
k2  (3.27  in  place  of  3.23  with  -F3  on  B(u)  € (0,4.42),  (5.84,6.43). 
Additional  improvement  is,  of  course,  still  possible  because  -F  does  not  lie 
in  B (u>)  , for  all  w . 


3.4  Case  Four  - Medium  and  High  Frequency  Modification. 

In  this  case,  P has  three  rhp  poles  and  two  rhp  zeros  with  coefficients 
such  that  P^  is  as  shown  in  Figure  9a,  i.e.,  in  the  first  quadrant  for  some 
(0,01^)  , in  the  fourth  for  (u^,u2)  and  back  in  the  first  for  (u 2,»)  , with 
|P(jco2)|  > |P(0)|  > IPfjoijH  and  P<«)  - 0 , e.g., 

2 


k P, 


-k(s  - .8) 


(s  - .2)  (s-6)  (s-10) 


, k 6 [l,k_]  . since  a constant  -F  is  capable 


, . „ . ..  , ....  _ K(s  + .2) (s+6) (s+10) 

of  handling  some  range  of  k , it  is  reasonable  to  try  -F.  - = , 

(s  + .8) 

choosing  K to  maximize  k2  ; giving  k2  « 4.36,  K ■ 351.6  . The  root  loci  of 


(1+P/F)  are  shown  in  Figure  9b  (qualitative) . Note  that  -F^  does  not  lie 
on  B(u>)  in  (4.0,14.1)  and  (37. 4,*)  and  the  relations  between  these  points 
in  Figures  9a, b.  Hence,  0_p  can  be  changed  in  these  intervals,  to  permit 
larger  k2  . 

How  should  0_p  be  modified  in  these  two  intervals?  In  this  case 

£ .■  -1.5  is  required,  so  there  must  exist  u>  ,u>  > w ,3  F(0)  < P (0)  , 

r x y X 1 

(f< j«*»x)  i > k2lp1<3“x)l  * < 'Pl(jV'  ' F<")  ” 0 with 

9_p  * 0O  at  u>  , to  . (For  F.  , to  » .75  , to  - 4 . There  is  another 
* r x y xx  y 

to  , to  pair  at  (14.1,  37.4)  because  F(“)  » • but  they  need  not  be 
x y 

considered.)  In  Figure  9a,  it  is  seen  that  | F( jtox)/F(0) ) is  the  constraining 

factor,  so  the  freedom  in  0_p  in  the  two  intervals  should  be  used  to  maximize 

this  factor.  From  the  sign  of  u(to)  in  (8),  this  means  decreasing  0 in 

these  two  intervals,  because  tox  < 4 . In  the  first  interval,  this  can  be  done 

2 2 

by  letting  F2  - K^ts  +B)/(s  +A)  , A < B and  choosing  K2  , B , A to 

maximize  k2  . When  this  is  done  an  improvement  is  indeed  found,  but  it  is 

very  small  (4.375)  and  the  region  with  -F  on  B(<o)  increases  to  ~ (0,8-9) 

and  (27.8,58.1)  . In  the  second  interval  the  modification  to  F would  logically 

2 

consist  of  a factor  (s  +B)  in  the  denominator.  This  leads  to  -F3  on  B(w) 
in  (0,10.6)  and  k2  ■ 4.379  , with  the  root  loci  in  Figure  9c,  for  this 
modification  only.  One  could  use  both  modifications  but  in  view  of  the  very 
small  k2  increases  for  F2  , 


F3  the  effort  is  not  worthwhile. 
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4.  CASE  CONSTANT  F NOT  SATISFACTORY 

In  all  previous  problem  classes  there  existed  a finite  k range  for 

which  -F  a constant  was  satisfactory.  This  was  due  to  -P  curving 

sufficiently  such  that  it  encircled  negatively  some  point  on  the  real  axis 

N/2  times,  with  N the  number  of  rhp  poles  of  P . Our  technique  was  based 

on  this  assumption,  as  one  could  then  start  with  0_p  » 0p  giving  -F^  the 

jto  axis  image  of  P . The  resulting  F^  could  be  improved  if  there  existed 

any  u>  interval  in  which  -F  did  not  lie  on  the  B(w)  . This  is  not  always 

the  case,  e.g.,  P = (s-4^~  has  6p  < 0 * <i>  € [0,«°)  , giving 

£ • - .5  at  most  for  -F  a constant.  Use  of  0_p  * 0p  no  longer  gives  a 

satisfactory  design  for  any  finite  k range,  no  matter  how  small.  It  is  shown 

below  how  any  such  P cm  be  modified  to  a new  P which  does  have  the 

desired  B.  = - N/2  for  -F  a constant.  One  cm  then  find  F for  this  new 

i| 1 sup 

P as  before.  Furthermore,  optimization  of  the  new  P is  also  possible,  as 
illustrated  later.  The  modification  of  p involves  insertion  of  additional 
rhp  ze  os  and/or  poles. 

n n n n 

Let  P ■ t p u * - ~ , be  the  modified  P and  T •=  r—  - -r—z — — =* 

d old  d d 1+P  d d + n n 

o 00 

n no 

— - — where  the  n^^  , d^  are  polynomials  with  real  coefficients.  The  desired 

N/2  encirclement  condition  for  P is  equivalent  to  D having  no  rhp  zeros. 

If  D has  only  interior  thp  zeros  then  continuity  of  the  zeros  of  D with 
respect  to  its  parameters,  guarMtees  the  existence  of  a finite  range  of  each 


for  which  D has  only  Ihp  zeros.  At  each  rhp  zero  of  dQ  (denoted  by  zi  , 


>m)  , T - 1 


m . 

and  n(s)  « £ A s1  is 


chosen  to  achieve  this  as  follows. 


i - 0 
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Let  D(s)  be  a real  Hurwitz  polynomial  of  degree  ? , larger  of 

[degree  of  d (s)  , or  of  n n (s)  ] , none  of  whose  zeros  are  those  of 
o o 

n (s)  but  is  otherwise  arbitrary.  Setting  n(z.)n  (z.)  - D(z.)  gives  m+1 

1 O 1 1 

linear  equations  in  the  m+1  coefficients.  As  the  resulting  matrix  in 

the  Vandermondian , there  is  a solution,  giving  n(s)  . Then  d(s)  is  obtained 

from  d d I ■ D-n  n , in  which  D , n , d , n are  known.  In  the  above 
o o o O 

the  rhp  zeros  of  dQ  were  assumed  distinct.  One  can  extend  the  proof  to  the 
general  case  but  it  is  simpler  to  suggest  that  one  can  always  perturb  the 
multiple  root  case  to  distinct  roots  with  negligible  effect  on  the  problem. 

The  above  procedure  is  certainly  practical  but  one  may  decrease  the  degrees 
of  n(s)  , d(s)  by  a suitable  rather  than  arbitrary  choice  of  D(s)  , as 
shown  in  a later  example. 

The  above  constitutes  proof  of  the  following  theorem.  Its  generalization 
to  non-real  P , G and  to  an  arbitrary  half  plane  is  obvious.  It  is  equally 
obvious  that  it  may  be  extended  to  apply  to  nonzero  independent  variations  in 
the  poles  and  zeros  of  P , as  well  as  to  its  gain  factor  (cf  2.  ). 

no<s> 

Theorem:  Given  P(s)  ■ k ^ , with  nQ  , dQ  real  polynomials  with 

o 

arbitrary  zeros.  It  is  always  possible  to  find  a real  rational  function  G(s)  , 
such  that  l+P(s)G(s)  has  all  its  zeros  in  the  interior  of  the  ihp  for  a 
nonzero  range  of  k . 

Example.  P - , for  which  constant  P is  no  solution  as  P(ju>) 

is  in  the  fourth  quadrant  V u € (0,*»)  . Using  the  previous  notation,  try 

first  T - (,~2)  and  see  if  a,B  >0  can  be  chosen  3 T(3)  - T(4)  - 1 . 

8 +As+B 
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This  is  not  possible.  It  is  also  found  that  in  T ■ — |M3.  AUS~2J — 

(s  +As+B) (s+C) 

it  is  impossible  to  find  A,B,C  >0  3 T(3)  • T(4)  ■ 1 . However,  it  can 

be  done  with  T - k(s~1^<3~2)  (-8:fp)  > e.g.,  A - B - 1 , p - -31/6  , K = -3  . 

(s  +AS+B) 


Thus  for  P 


new 


P*  - 


k(s-l)  (s-2)  (s-^) 
(s  -|)  (s-3)  (s-4) 


there  exists  a constant  P , 3 f+p* 


is  Hurwitz  for  a finite  range  of  k . 


Optimization  is  achieved  by  first  letting  9_p  ■ 9p<l  giving 


Kx(s  +|)  (s+3)  (s+4) 
(s+1)  (s+2)  (s  +^-) 


so  that  the  root  loci  of  1 - P*/F  are  as  shown  in 


Figure  10  with  k^  - km^n  ■ the  larger  of  k^ , k^2  and  k2  the  smaller 

of  (l,k21>  . An  optimum  design  would  be  achieved  if  k^  - k12  and  k21  - l 
(giving  roots  at  ± j«»)  , for  then  -F(jw)  lies  on  B(u)  V u £ |0,»)  . 


He  therefore  .rite  P*  - 1*:fl  . F - !**"  j**4’  end 

(s-p) (s-3) (s-4)  (s+1) (s+2) (s+z) 

seek  z , p to  achieve  the  two  conditions  k^  ■ k^2  , k21  * 1 . The 

result  is  z - 5.275  , p - .832  , giving  ■ 1.117  (in  place  of 

1.068  for  z * 31/6  , p ■ 5/6  ) and  -F  on  B(w)  V u € t0,«»)  . 


>•! 

% 
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5.  CONCLUSIONS 

A design  philosophy  and  methodology  has  been  presented  which  is 
useful  for  eliminating  or  reducing  the  nmp  property  for  structures  of  the 
form  of  Figure  1.  One  first  deals  with  a P function,  such  that  F a 
constant  guarantees  the  Hurwitz  character  of  the  numerator  of  1+P/F  , over 
some  nonzero  range  of  the  gain  factor  k of  P . It  was  shown  in  4.  that 
the  original  P can  always  be  modified  to  have  this  property.  Then 

“ ®p  gives  a first  attenpt  F^  , which  gives  a larger  tolerable  range 
of  k than  F a constant.  If  the  resulting  -F  does  not  lie  on  b(u)  V u>  , 
then  F can  be  improved . The  means  of  improvement  lies  in  deciding  whether 
additional  phase  lag  or  lead  of  F is  desirable  in  those  intervals  where 
~Fj  does  not  lie  on  B(u)  . The  Bode  integrals  are  used  in  making  this 
decision  and  the  implementation  is  by  means  of  poles  and  zeros  on  the  jw 
axis,  such  that  0p  is  unchanged  in  the  other  intervals. 

However,  the  optimization  of  F was  always  over  the  irps  class,  and 
it  is  conceivable  that  non-mps  F may  be  superior.  One  way  to  handle  this 
possibility  in  the  context  of  this  paper,  is  to  imbed  P in  a larger  nmp 
unstable  class  with  more  rhp  poles  or/and  zeros  than  the  original  and  apply  the 
approach  of  this  paper.  'Riese  additional  rhp  critical  points  are  of  course 
supplied  by  F . Thus,  direct  optimization  on  the  most  general  class  of  F 
is  not  presented  in  this  paper. 
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Uncertainty  in  the  poles  and  zeros  of  P is  easiest  to  handle  when  these 
uncertainties  are  independent.  For  then,  there  exist  extreme  pole-zero  com- 
binations which  limit  the  maximum  range  of  k and  give  the  resulting  F . The 
latter  is  then  satisfactory  for  all  the  other  pole-zero  possibilities.  These 
extreme  pole-zero  combinations  may  be  such  that  no  F is  at  all  available. 

The  problem  of  finding  Fgup  can  be  much  more  difficult  if  the  uncertainties 
are  related. 
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